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Abstract. We investigate the spin dynamics of excitons interacting with 
acoustic phonons in quantum wells, quantum wires and quantum disks by 
employing a multiband model based on the 4x4 Luttinger Hamiltonian. 
We also use the Bir-Pikus Hamiltonian to model the coupling of excitons to 
both longitudinal acoustic phonons and transverse acoustic phonons, thereby 
providing us with a realistic framework in which to determine details of the spin 
dynamics of excitons. We use a fractional dimensional formulation to model 
the excitonic wavcfunctions and we demonstrate explicitly the decrease of spin 
relaxation time with dimensionality. Our numerical results are consistent with 
experimental results of spin relaxation times for various configurations of the 
GaAs/Alo.3Gao.7As material system. We find that longitudinal and transverse 
acoustic phonons are equally significant in processes of exciton spin relaxations 
involving acoustic phonons. 



PACS numbers: 73.21.Fg, 72.25.Fe, 78.67.De, 73.63.Hs 

1. Introduction 

Excitons in low dimensional systems, formed due to the binding of electrons and holes, 
have well known features such as enhanced binding energies and oscillator strengths 
. In recent years there have been increased studies of excitonic processes involving 
phonons within the restricted dimensional exciton energy band and exciton spin 

relaxation jl] , formation and phonon-assisted thermalization El Ej have attracted 
great interest. This is largely due to the significant advancement in spectroscopic 
techniques involving continuous-wave and time-resolved photoluminesccncc excitation 
measurements 8 s which provide results of high temporal and spatial resolution. In 
particular, the use of polarization dependent photoluminescence 110] involving 
circularly polarized light has extended the possible application of excitonic systems 
to spintronics |1 II IT2] in which the spin degree of freedom is exploited in electronic 
devices. 

Several processes involving the Elliot- Yafet mechanism, the Bir-Aroniv-Pikus 
mechanism and the Dyakonov-Perel mechanism are known to cause spin relaxation 
of carriers in solids ^3 . We are interested in the first of two other possible 
mechanisms which can cause exciton scattering involving spin-flip processes, firstly, 
holes interacting with acoustic phonons ^JEDEH anci , secondly, long-range electron- 
hole exchange interactions which couple the spins of electrons and holes |17| . In narrow 
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wells both mechanisms can become comparable, the dominant process depending on 
the electron density, temperature and quantum well width. 

Spin dynamics of excitons are identified through measurements of circular 
polarization [5|-|1U| in which, initially, spin polarized electrons and holes are formed 
followed by both spin and energy relaxation of carriers and, finally, there is 
recombination of carriers resulting in the measurement of a circularly polarized signal. 
It is not clear, however, whether the electron and hole undergo spin and energy 
relaxation as unlinked carriers and then become coupled (forming an exciton) after 
which they recombine, or whether the carriers become coupled first, then undergo spin 
and energy relaxation before the recombination process. Whether there is concurrence 
of both spin and energy relaxation in electron-hole systems also needs careful study. 
There have been suggestions ^H] that spin relaxation occurs after energy relaxation 
into the excitonic ground state, while others ^Hj have suggested that momentum 
relaxation occurs simultaneously with spin relaxation in excitonic systems. This 
interplay of energy and spin relaxation processes needs further study in order to gain a 
thorough understanding of phonon related dynamics in low dimensional semiconductor 
systems. Such details are also vital for the accurate interpretation of experimental 
results. 

In a simple exciton model the fourfold degeneracy at the top of the valence band, 
which introduces some nontrivial features in the excitonic spectrum, is neglected. 
The complicated nature of the multiband Wannier exciton model was noted by 
Dresselhaus who demonstrated the absence of a well-defined exciton centre-of- 
mass transformation. Subsequently Baldereschi and Lipari [23 |^ studied exciton 
dispersion relations for a variety of semiconductor systems using a generalized centre- 
of-mass transformation. In the Baldereschi-Lipari approach, the "d-like" term is 
treated as a perturbation on the "s-like" term of the exciton Hamiltonian which has 
a hydrogen-like spectrum. Using this approach, Kane [22 obtained useful expressions 
for the energy shift of the bulk Is exciton from the conduction band minimum. In 
recent years, the multiband model of the exciton interacting with phonons in bulk 
semiconductor systems, investigated in [231 has been extended to the quasi-two 
dimensional model of the exciton [1511251 05) . These papers show that the description 
of low-dimensional excitons within a multiband model is a complex problem that 
requires extensive computational efforts. Here we use a fractional dimensional 
formalism to extend the study of spin dynamics of excitons to quantum wires and 
quantum disks and to emphasize the significance of the effect of dimensionality in the 
spin dynamics of excitons. 

The paper is organized as follows. In Section [2] we present the theoretical 
formulation of the quasi-two dimensional exciton model that is used to obtain suitable 
dispersion relations as well as an explicit form of the exciton wavevector. In Sectional 
we give details of the fractional dimensional approach and present results specific to 
the quantum wire system. In Section0|we use the Bir-Pikus Hamiltonian to model the 
coupling of excitons to acoustic phonons and obtain expressions of matrix elements 
corresponding to longitudinal acoustic (LA) and transverse acoustic (TA) phonon 
modes. In Section we extend the theory developed in Sections El and 01 to obtain 
expressions of spin transitions corresponding to the decay of heavy- and light-hole 
excitons and intraband as well as interband scattering in constrained systems. We 
report and discuss the results of salient features of excitonic processes mediated by 
LA and TA phonons in Section H3 and summarize our conclusions in Section 
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2. Theoretical Formulation 

We investigate properties of Wannier excitons in quantum wells fabricated using direct- 
gap cubic semiconductors, using as an example GaAs which is representative of a 
large class of III-V and II- VI zinc blende semiconductors. The conduction band edge 
in such systems involves an s— type spin degenerate carrier state while the valence 
band edge involves a p— symmetry state with a sixfold degeneracy [21 E3 - When spin- 
orbit interactions are included the sixfold degenerate level F 5 splits into a fourfold 
degenerate state r 8 (with J — |) and a twofold degenerate state T 7 (J = i). We 
neglect the twofold degenerate state, also known as the splitoff band, in this paper. 
When analysing these spin states, we use a basis that is oriented relative to the 
quantization axis, where the Z-axis is chosen in the growth direction of the quantum 
well structure. The choice of a preferential axis of quantization for the hole spin 
subsequently leads to notable anisotropic properties of the exciton as we show later 
in this work. 

The vectors in the angular momentum basis generated by J = (J x ,Jy,Jz) 
projected on the quantization axis are labelled by the magnetic quantum numbers 
m = ±| for the heavy hole (HH) states and m = ±i for the light- hole (LH) states, for 
the degenerate level r§ (see Chapter 3 in 0). The broken translational symmetry in 
the growth direction of quantum wells lifts the degeneracy of the heavy and light-hole 
valence bands at the T-point, giving rise to a distinct splitting at the two-dimensional 
in-plane hole wavevector = 0. This is the basis of the formation of heavy- and light- 
hole excitons in quantum wells. At k^ ^ 0, the LH and HH states mix and hybridize 
giving rise to strong nonparabolicity [211123 m the valence band spectrum. Due to the 
distortion of the valence band structure, we classify the heavy- and light-hole exciton 
bands by their basis functions at — 0. Unlike the split-off band which arises from 
the spin-orbit coupling, the splitting between the heavy-hole and light-hole states is 
a consequence of strain due to confinement in the growth direction. Throughout this 
paper k c (feh) marked with an arrow denotes the in-plane electron (hole) wavevector, 
while k c (kh) denotes the three-dimensional electron (hole) wavevector. 

In order to define operators associated with the spins of charge carriers in quantum 
wells we use the standard rotation matrix D(a, (3, 7), which is a function of the Euler 
angles (where < a < 2tt, 0^/3^ 7r, ^ 7 < 2tt), and is given by [3*Tll3^] : 

D(a,j3,^) — exp(— ia J z ) exp(— i(3J y ) exp(— i7 J z ). 

The angles a = <fi and (3 = 9 are in fact the polar angles of the spin state J(6, </>) with 
respect to the frame in which Z is the quantization axis. A rotated spin state can be 
represented as a combination of the 2 J + 1 component states using standard rotation 
matrices: 

D(a,p n )\Jm) = ^2(Jm'\D(a,l3 n )\Jm) \ Jm!). (1) 

The matrix elements D J m , m {a, (3, 7) = (Jm'\D(a, (3, 7)! Jm), where — J ^ m,m' ^ J, 
may be written out explicitly for any J, however we require only the explicit form and 
transformation properties for angular momentum values J — \ and J = | which are 
applicable to the electron and hole. Details are given by Brink and Satchler [32]. As 
the basis vectors of the representation in equation Q are chosen as the eigenfunctions 
of J z , we use the simplification |D^, m (a, /?, 7)! = |d„' m (/3)| where d J is the reduced 
rotation matrix. A variety of relations among the matrix elements d^^fi) is given 
in Chapter 3 of the monograph by Biedenharn and Louck |33j . 
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2.1. Electron and hole operators with arbitrary spin orientation. 

By using properties of the reduced rotation matrices we are able to write the quasi- 
two dimensional electron and heavy- (or light-hole) creation operators as linear 
combinations of operators corresponding to the spin projections on the quantization 
axis. The creation operators of the electron at (5), heavy hole ct (77) and light-hole 

ct (77) spinor systems thus appear as follows: 

&z 

where t z = ±5 and a z — ±§(±|) are the respective electron and heavy (light-) hole 
spin projection eigenvalues of J z . The indices <; = ± and 77 = ± denote the two 
possible degenerate states of the electron (at fc c = 0) and heavy and light-holes (at 
fch = 0) respectively, and 9 e (9h) denotes the azimuthal angle of spin states of the 
electron (hole) operators with the in-plane wave vector fc (fch)- For simplicity, we use 
the same notation for the heavy- and light-hole wavevectors and the corresponding 
azimuthal angles. 

The dependence of the operators at (<^) , ct (77) and ct (77) ) on their respective 

azimuthal angles yields the expected relations at selected angles. For instance, 9 C = 
yields the relations: at (+) = at (i) and at (— ) = at {—\). Likewise 9 C — ir yields 

the relations: at (+) = —at (— i) and at (— ) = —at (i). Similar relations can be 

rZ c /Cq ftj<2 

obtained for the two types of hole spinor systems. 



2.2. Exciton Wavefunction and Spin 

The field operator ^ of an exciton, with a centre of mass that moves freely, is 
composed of electron and hole operators with spins aligned along the quantization 
axis, and has the expression: 



*t(re,rh)= E 



4» e ifc -- ^< : (r c ) a\{r z )a^{r z ) (3) 

Tz 



where r c = (x c ,y c ,z c ) and rh = (x^, y^, Zh) denote the space coordinates of the 
electron and hole, respectively, and f c (or r\) marked with an arrow denotes the 
in-plane coordinates of the electron (or hole). u^(r e ) (wl^(rh)) is the periodic 
component of the bulk Bloch function of the electron (hole) which depends on the 
projection t z (<j z ) of the spin momentum onto the positive Z-axis. As usual we 
assume that the Bloch functions are identical in both well and barrier materials. The 
conduction and valence bands in semiconductors with zinc blende structures are made 
up from s- and p- type atomic functions respectively. Denoting now by | U) the orbital 
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component of the s-type Bloch functions, we write the electron Bloch wavefunctions 
as: 

«kf(rc) =|§,£> = |tft>, 

< /2 (r c ) = |i,-i) = |E/j), (4) 

where as usual denote spin-up and spin-down states. Likewise, by denoting the 
Bloch functions with the symmetry of p x , p y , and p z orbitals as \X), \Y), and \Z), 
respectively, we may write the hole Bloch wavefunctions as linear combinations of 
\(X + iY)i), |(X + iY)T), \(X-iY)l), |(X-iy)T), \Z{) and \Z]). 

It is well known that the conservation of the in-plane wavevectors fch and fc e in 
quantum wells allows the lateral motion of the charge carriers to be decoupled from 
their motion along the Z-axis. This gives rise to subbands of quantized envelope 
functions, £r (z c ) and £^(zh), in the conduction and valence bands respectively. 
The electron subband function which has even or odd parity under z e inversion is 
determined by using the eigenfunction of the BenDaniel-Duke model |34) in which the 
potential is approximated using a finite square well. The boundary conditions at the 
interfaces are determined by using the wavefunction and current continuity relations 
and 

E «"■"*• i (5) 

where fc c = (fc c , <j) c ). 

Calculation of the hole envelope function £,°*{zh) is not as straight forward, as 

its dependence on fch and spin components o~ z gives rise to a complicated warping of 
the various valence subbands |30j . The axial approximation, in which a cylindrical 
symmetry around the crystal growth direction is assumed (following |15p. simplifies 
the hole subband states which may be expressed as a superposition of four states of 
the valence band top with spin projections o~ z = §, 5, — 51 ~§ : 

= £^"*fc|(^ (6) 

where fch = (fch,</>h). Further details of the evaluation of the hole wavefunctions, 
characterized by the Luttinger parameters 71,72 and 73, are given in the Appendix. 

An exciton carrying a centre-of-mass momentum K is also characterized by its 
spin S and its projection S z on the quantization axis. A given exciton spin state 
depends on the spin-spin correlations which exist between the electron and hole spinor 
systems. We consider only the Is exciton so that the exciton spin S is determined 
by the group theoretical rules of angular momentum vector addition. Thus the four- 
dimensional state space of the exciton system spanned by the states \S) is characterized 
by a coupled system comprising the two-dimensional electron (|r z )) and hole (|cr 2 )) 
spin states. In this paper we choose the exciton basis states with respect to the same 
axis of quantization as that of the hole and electron. 

We define a generalized exciton as one that involves the coupling of an electron 
with wavevector fc e and spin orientation angle 9 C to a heavy- or light-hole with 
wavevector fch and spin orientation angle 9^. At low densities the exciton creation and 
annihilation operators satisfy boson commutation relations, but at higher densities we 
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write the exciton eigenvector \K\S) explicitly as a linear combination of all products 
of possible pairs of electron and hole eigenvectors: 

\K;S) (7) 



^2 \Fu{k c ,kh-S,V) 



k e ,K;<;,r)) + F h (k c ,k h ;<;,'q) 



kcjk^: 77 



where Fn(k e , fch! ?, 77) (Fh(k e , kh,; <;, 77)) represents the coefficient function of coupling 
between the electron and heavy-hole (light-hole). We have omitted the azimuthal 
angles 9 C and 9h of the spin states in equation iJJ}, for simplicity in notation. The 
eigenvectors \k e , k^; <;, r])yi and |fc e , fch; v)l are constructed from linear combinations 
of the electron and hole creation operators of the respective forms: at (5) ct (77) |0) 

and at (<^) ct (77) 1 0) , where |0) is the vacuum state. We have considered only the 

lowest heavy-hole and light-hole exciton states and therefore neglected energy band 
indices in equation Q. The exciton wavevector in equation J7J contains spin mixing 
effects between heavy and light- hole exciton states, which forms the basis by which 
acoustic phonons contribute to spin relaxation processes and will be further studied 
in Section 01 



2.3. Exciton spin states at 9 C = 9^ = 

At 9 C = 6>h = the electron and hole spins are oriented in the direction of the 
quantization axis and specific forms for the electron and hole creation operators are 
given by the following, firstly for heavy-hole excitons: 



\k e , fch! 5) H 



E 

T»=±- 



and, secondly, for light- hole excitons: 
|fe e ,fch;5) L = 



E 



,2 Tz ' 2 ( 



±T • iff 



55 2 )4 o (r.) 4 h H (ff,)|0) (8) 



SS z )4 c (r z )4 hL (ff z )|0), (9) 



tr,=±A 



where |0) is the electronic vacuum state of the system representing completely filled 
valence band and empty conduction bands, and the factors (^t z ; §ff z |55 z ) and 
(\t z ') h°~z\'SS z ) are Clebsch-Gordon coefficients. By using equations (JSJl and © wc 
obtain the following exciton spin states (for 5 = and 5 = 1) at 9 C = 9^ = 0: 

|fc e , 4; S = 0) L = ^ [at, (J) 4 J- 1) - a^-i) c^i)] |0>, (10) 
l^4;S= 1) L = £ [atj-i) 4 h)L (|) + at e (i) ^ J-i)] |0) 

+^[4 e (-i)i )L (-l)+4 e (i)i )L (i)]i Q )' (n) 

1^5 5 = Dh = 73 + 4 C (5) 4 h , H (-D] l°>" ( 12 ) 

The 5 = spin state of the exciton arises from the coupling of the electron to the 
light-hole while the 5 = 1 spin state results from the combination of three symmetrical 
spin functions of the electron and light-hole spin states with the spin components 
S z = +1,0, —1. The 5 = 1 exciton can also be formed by coupling the components 
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S z = ±1 between electron and heavy-hole spin states. Spin states S = 2 of the heavy- 
hole exciton cannot couple to photons in the light field and are therefore known as 
optically inactive dark excitons. Hence, the light-hole and heavy-hole excitons are 
characterized by the spin projections S z = 0, ±1 and S z = ±1, ±2 of the total angular 
momentum S = 1,2 respectively. The short range exchange interaction (given below 
in equation (jl8(l ^ further splits the ground states of both LH and HH excitons into 
doublet states, an effect also known as singlet-triplet splitting (see jSS] for a recent 
review). Hence, in low dimensional systems spin mixing effects between heavy- and 
light-hole states result in a complicated energy spectrum and the spin states for S = 
or S = 1 are no longer eigenfunctions of the corresponding Hamiltonian. In this paper 
we neglect heavy- hole and light-hole mixing effects due to exchange interactions. We 
also assume that the mutual conversion between the heavy-hole and light-hole exciton 
states arises from the phonon scattering processes as discussed in detail below in S^l 



2.4- Exciton eigenvalue problem 



The exciton eigenvalue E CX (K; S) corresponding to the eigenvector \K ; S) in equation 
Q is found by solving the Schrodinger equation: 

H loc + Hg t + H^ h ) \K; S) = E ex (K; S)\K; S), 



(13) 



where the Hamiltonian is written as a sum of four terms which we now describe in 
turn, flkin is the kinetic energy: 



kin 



E 



n 2 e 

En 

9 2m 


4 C W a fccW 






2m L 



(14) 



21,2 



where m c is the effective electron mass and win (tol) is the heavy-hole (light- hole) 
mass. 

I/ioc is associated with the localization energies of electron and hole, and is given 

by 



H\ oc — 

where 



2m r 



R> = y % (2 C ) (-ift^- J e £c (*,), (16) 

with a similar expression for (jrL), and L w denotes the thickness of the quantum well. 

HP t denotes the Hamiltonian due to the Coulomb interaction which scatters 
an electron in the initial state with momentum k and spin ^ to a final state with 
momentum k + q and spin <;' . Likewise a hole in its initial state with momentum h! 
and spin r\ is scattered to a final state with momentum k' — q and spin rf , and hence: 



H. 



\ E (k + q,^'; k'-q,rj' U(\r e - r h |) fc',77; fc,<r 



(17) 



k,k f ,q,X 
c n c> n' 



Xa l. + q^ 4>_* M) C P.X(V) Ojf(?) 



k'-q,X 
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where 



U(\r e -r h \) = - r , 

e|r c - r h | 

in which e is the relative dielectric constant and X denotes either the heavy-hole 
operator H or the light- hole operator L. In the two-particle interaction matrix element 



{K + q,^; K' 



[/(|r e -r h |) K',r); K,s) , 



the states to the right of the scattering potential U(\r e — rh|) represent the initial 
states while those to the left represent the final scattered states. 

H^ h denotes the Hamiltonian due to the exchange interaction which scatters an 
electron with spin ^ and momentum k! in the conduction band to a hole state carrying 
momentum k' — q and spin r\ in the valence band. At the same time, a hole state 
with momentum k and r\ in the valence band is scattered to an electron state with 
momentum k + q and spin in the conduction band: 



^fnt Ch = \ E - V> k + q^ U(\r c -r h \)k,rj; k' , C 



(18) 



fc,fc',q,X 



xc L„- y W a L?M c k.x(v) H'^)- 



By using equations I|2I5I14I17I18|I and by setting 6 C = 6^ = we obtain, after 
some manipulation, the following two coupled equations for the unknown coefficient 
functions Fn(k c , k^\ rf) and F\^(k e , fch; S", rf) which are defined by equation Q: 







E e (k c ) + E H (k h ) - E*(K) F H (K, k; <;, v) 



(19) 



+ E / ( V c& k'; H)-(l- S) V cxch (k, k'; ff)) 



x 
and 
= 



F u (K,k';,,r ] )d 1/ 1 2 T (<f> e ) d 6 ' 3 \ (fa) + F h (K, k'; s, rj) d 1 " (fa) 

S o T z r l o ° z 2 Tz r l 2 ° z 



j3/2 



1/2 



i 3/2 



E c (k c ) + E L (k h ) - E^(K) F L (K, k; ? , V ) 



(20) 



F K (R,i?-,S,Tl)<t{ 2 T (0h) + FL(^,fc';^ ? ?) dT r (^Kl^ (</>h) 

^2 'z 'I 2 z ^ 2 z I 2 z 

where E c (k e ) and £Jji(^h) (^L(fch)) denote the in-plane energies of the electron and 
heavy-hole (light-hole) respectively, and E^ X (K) (E^ X (K)) denotes the heavy-hole 
(light-hole) exciton energy in the plane of the quantum well. The angles <p c and 
0h are defined by equations JSJ) and ©. 

The integral forms of Vc(k, k' ; X) and T4 xc h(fc, k! \ X) which are the in-plane 
Fourier transforms of the Coulomb potential, with X = H (a z ,a' z = ±|) or X = L 
(a z , a' z = ±±), are given by: 

V c (k,k';X) = ? 1 =r J2 e 1 ^^"^^) (21) 



j3/2 



1/2 



j3/2 



27T6 \k - k' 
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V exch (k,k';X) = J2 e^^~<^ (22) 

Equations (|19fl and l|2(J|l give rise to strong nonparabolic exciton centre-of- 
mass dispersions and require a large basis function set and sophisticated numerical 
techniques in order to be solved, even for the particular case 9 e = 9^ = 0. As a 
result of the nonparabolicity there is no unique definition for the mass of the exciton, 
which has resulted in a range of mass values in the literature |36| . The hole spin 
components also become mixed and some authors |15| have found it convenient to use 
the parity-flip mechanism to incorporate changes in parity in the valence band states. 

The centre-of-mass momentum wavevector K and the relative wavevector k are 
defined in terms of the electron and hole wavevectors (fc c and fch respectively) by: 

K = k G k\ L , 

k = (1-1?) ke + ti 4, (23) 

where the parameter $ is generally used in a scalar form (5l l3fill3"?| which conveniently 
allows the decoupling of the exciton centre-of-mass motion from its relative motion. In 
more rigorous calculations involving the heavy-hole and light-hole dispersion relations, 
d has been used as a tensor |2S1 Recenty, Siarkos et al 25, used a variational 
approach involving minimization of the exciton energy and obtained optimized values 
of i9 in the range ^ i9 ^ 0.5 at \K\ ^ 0.4 nm -1 . For our work, we select suitable 
values for $ for the heavy- and light-hole Is exciton by extending the expression 
derived by Kane [221 for the energy shift of the bulk Is exciton from the conduction 
band minimum to a low dimensional K— space: 

— = — - — (E) 2 2 

M a M 3 331 )S 72 

1 2s 2 2 

IF =—72 (24) 
M c mo 

where the + sign and X = H correspond to the heavy-hole exciton while the — sign 
and X = L correspond to the light-hole exciton. R y is the effective Rydberg constant, 
m is the free electron mass and 

-i ™>c 1 
M =j 1 +m , s = ; — r. 

The functions g\{E) and gs{E) given by Kane [221 are assumed to apply in low 
dimensional systems, and the constants 71,72,73 are discussed in the Appendix. By 
using E = —R y and the approximation 73 w 72 we obtain: 

q TO « . 2 L 1 5 ^ 2 
a c = v = ——: h 2s 72 m ±1 —r 

7l + TO C L 71 + m c 

at = 1 - ti. (25) 
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With the value m c = 0.065 m for the effective electron mass and 71,72 as given in 
table^for the GaAs/Al^Gai-^As material system, we obtain a c — 0.377 for the heavy- 
hole exciton and a c = 0.116 for the light-hole exciton. These values are consistent 
with those obtained in reference |25) . 

In order to obtain an explicit form for the exciton wavevector, we neglect the 
relatively weak exchange interaction term T4 X ch(&, k' , X) involving an overlap of the 
conduction and valence bands. We now write equation J2} in the form: 

\K;S) = |A J ;5) H + |/?;5) L (26) 

and use equations (|19fl and 120|) to obtain (for the heavy-hole exciton): 

\K- 5) H = _ £ E W <f W 6 KA-Jk *is(a c k h + a h k c ) 

xfcfe)fhfe) 4 o (r,) 4 h (a z )|0), (27) 
and (for the light-hole exciton): 

\K;S) L = ^ £ * ls(ac fc h + ah fc c) 

*«=±5 

xFcfe) Fh(fc) 4 (t,) 4 (o-*)|0) (28) 

where is the angle of scattering and /i denotes the polarization state of the exciton. 
We have also dropped the X = H(L) label from the hole creation operator as we have 
restricted the azimuthal angles to 6 C = 6h — 0, and therefore use the same notation 
ct(c z ) for the heavy-hole and light-hole creation operator. 
The form factors F c (q z ) and Fh(q z ) are given by: 



F c (q z ) = j dz c e^ z c ^JZe) 

F h (q z )= J dz h e^^|^ h| (z h )| 2 - (29) 

The exciton wavevectors in equations l|26Jl - (|28(l generalize the form used in earlier 
work |H] in which we assumed a simple two-band model without taking into account 
the fourfold degeneracies of the upper valence energy bands. \E , i s (a e ^h + Qh^c) is the 
wavefunction of a hydrogen type system in momentum space, which depends on the 
relative electron-hole separation. This function determines the intrinsic properties of 
the exciton and forms the basis of the fractional dimensional approach detailed in the 
next section. 



3. The fractional dimensional formalism 

In the fractional dimensional approach the confined "exciton and quantum well" , 
"exciton and quantum wire" or "exciton and disk" system is modelled using 
conventional Hilbert spaces which carry a fractional dimensional parameter, and 
in which the composite exciton system behaves in an unconfined manner, and so 
the fractional dimension is related to the degree of confinement of the physical 
system. Such an approach introduces simplicity and utility and has been used in 
the investigation of several important processes in low dimensional systems |38 j -|42 | . 
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Very recent calculations in quantum wires |43j with moderate strength of coupling 
between the subbands have shown that the fractional dimensional approach provides 
strikingly good agreement with results obtained using models that require intensive 
computational efforts. In recent work |44| we have developed algebraic properties 
of various quantum mechanical systems using fractional dimensions. These studies 
demonstrate the useful and important role of the dimensionality (denoted here by a, 
following |48| ) which interpolates between zero in ideal quantum dots to three in an 
exact three dimensional system, such as an infinitely wide quantum well. 
The binding energy of a fractional dimensional exciton |45| is given by: 

E b = —. R \, , (30) 

where n = 1, 2, ... is the principal quantum number of the exciton internal state and 
R y is the effective Rydberg constant. The function ty\ s (k) appearing in equations (|27[1 
and H28|) . where k = a c kh + a^ke, has been obtained (Thilagam |38p in the form: 

vMfc) = (4vr)t-i(a- l)f+yr [l(a l)](a B )f (l + [±(c* - 1) \k\ a B ] J " (31) 

where a B is the three-dimensional Bohr radius of the exciton. \E'i s (A:) yields the 
expected forms in the exact three dimensional and two dimensional limits. The 
dimensional parameter a characterizes the degree of "compression" of the confined 
exciton and can be calculated in several ways (see for instance references |45] - 
471). Once a is determined, many useful properties (see |3H]-02], 02]) can be 
determined by substituting its value into the appropriate equation that corresponds 
to a physical quantity. Recently, Escorcia et al. @H] used variationally determined 
envelope functions to calculate the wavefunctions of the exciton relative motion in 
quantum wires and disks. Like earlier approaches [451 147) . the technique of evaluating 
a in reference |1SJ does not include valence band coupling effects between the heavy- 
and light-hole bands. For our work, we improve on earlier methods by utilizing 

realistic values of a that is evaluated by including valence band coupling effects through 
a 4 x 4 Baldereschi-Lipari Hamiltonian |2()l l2*T] which describes the relative motion of 
the exciton. This allows salient features due to the anisotropic exciton band model to 
be incorporated in a. 

We approximate the exciton ground-state wavefunction in quantum wires and 
quantum disks by equations i|26l28[l with appropriate electron and hole wavefunctions 
£|£ |(r e ) and |( r h) determined (following references [21 EH) by the quantization 

length R which is the radius of the quantum wire or quantum disk where the latter 
is also specified by a thickness L. The function \&i s (fc) is modelled using equation 
<|31[1 with predetermined values of a. In figure ^ wc compare the dimensionalities 
of the heavy- and light-hole excitons as functions of the quantum wire radius R in 
GaAs/AL E Gai_ ;E As using the parameters listed in tabled We have neglected the 
effective mass and dielectric mismatches between the host materials. 

At large R values the dimensionality of the light-hole exciton is noticeably smaller 
than that of the heavy-hole exciton. The dimensionalities reach minimum values at 
the strong confinement limit of R m w 3.5 x 10 -9 m. For R ^ R m the excitonic 
wavefunction spreads beyond the well material and leads to a rapid increase of the 
dimensionalities as shown in figure^ The effect of including the off-diagonal elements 
of the Luttinger Hamiltonian through the 4x4 Baldereschi-Lipari Hamiltonian |2()II21| 
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Figure 1. Dimensionalities of the heavy- and light-hole excitons as functions of 
the quantum wire radius R in GaAs/Al^Gai—zAs using the parameters in table 
m The dashed lines correspond to calculations in which off-diagonal elements of 
the 1 x 1 Baldereschi-Lipari Hamiltonian (references 1201 1211 1 are neglected. 
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can be seen in the difference between the soiid curves and dashed curves (in which 
the off-diagonal elements of the 4x4 Hamiltonian are neglected). This figure shows 
that the ideal lower limit of a = 1 is never attained. Moreover, because a ^ 2 
for R ^ 10 -8 m, dimensionalities in quantum wires are not restricted to the range 
1 ^ a ^ 2. These results also indicate that it is not appropriate to use equation (|30|l 
to estimate a in quantum wires. 



4. Exciton-phonon interaction 



Two main mechanisms are considered when analysing exciton-phonon interactions 
in semiconductors: the deformation potential and the piezoelectric interaction. We 
write the Hamiltonian H^_ ph describing the exciton-phonon interaction, due to 
the deformation potential, in terms of electron and hole creation and annihilation 
operators: 



H. 



DP 



-ph 



E 



M^(q,q z )ai + Jr z )a^T' z ) 



M^,q z )ct +? (a z )c^a' z ) 



k,q,q z 



-q z ) + b\(q,q z ) 



(32) 



where the phonon wave vector is given by q = (q x , q y ,q z ) = (q, q z ), and where b x (q, q z ) 
(respectively b\(q,q z ) ) denotes the creation (respectively annihilation) operator of a 
A— mode phonon. The A-mode is denoted LA for longitudinal acoustic phonons and 
TA for transverse acoustic phonons. The strain tensor due to the various A— mode 
phonons is written in terms of normal-mode coordinates as 



2pVu>qx 



+ b\(q,q z ) (eixqj + ejxq^e 1 ' 



(33) 



where p is the mass density of the material. The acoustic phonon energy spectrum 



is determined by ui 



\q\ for the longitudinal mode and to 



the transverse mode, with v LA and v TA denoting the corresponding sound velocities. 
The term e<> represents the unit vector of polarization of the A-phonon along the 
i-dircction. The longitudinal mode e LA and the two transverse modes e 
take the form ESI : 



I VI an °- ^TA2 




(34) 



where the TAl-mode (TA2-mode) corresponds to the transverse phonon polarized 
along the x (y)-direction. 

The matrix elements (q,q z ) and (q,q z ) in equation depend on the 
precise nature of the exciton-phonon interactions, such as the deformation potential 
interaction, the piezoelectric interaction, and the phonon mode. We have explicitly 
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included the spin indices to distinguish the hole states involved in the interaction. 
Due to the isotropic nature of the electron related interaction with acoustic phonons 
in cubic III-V semiconductors ESI, we obtain: 



2pVv 



LA 



Mli(q,q z )=0, (35) 

^ 2 

where S c is the deformation potential constant associated with the conduction band. 
In using equation (|3*5j) some papers (for example jKHJ) have neglected the deformation 
potential coupling to the TA phonon mode. The importance of including interactions 
due to TA phonons becomes clear in the context of hole-phonon interactions, which 
we consider next. 

The component of the hole related interaction in the exciton-phonon Hamiltonian 
is described by the Bir-Pikus Hamiltonian Hbp |511 15^] which takes the form 

/ F H I \ 

H* G / 

I* G —H 

\ I* -H* F J 

where the matrix elements are defined (following |51p in terms of deformation potential 
tensor components S^- and strain tensor components by: 

A = J2*t j e ij , (A = F,G,H,I) 

where 



H BP = 



(36) 





hi 










"11 


"F 
— "22 — 


a 


"33 


= b, 




n G 
"11 


"G 
— "22 — 


\{a + 2b) 


TjG 
"33 


= i(4a- 


-b), 


~H 
"13 


"H 

— "31 — 


— i c 


t;H 
"23 


— "32 — 


!"13' 


— / 
"11 


— / 

— "22 




^1 
"12 


—7 

— "21 — 


-ic, 



(37) 

with all other components 2^- = 0. The values a, b and c are evaluated using 
the matrix elements (X\E XX \X), (X\E yy \X) and (X\E xy \Y) + (X\E yx \Y) where the 
tensor operator £L- is obtained |52| by subtracting elements of the Hamiltonian of the 
unstrained crystal from those in a strained crystal. 

The effective deformation potential interaction can be obtained separately for 
TA and LA hole-phonon interactions via a suitable transformation of the deformation 
potential tensor matrix: 

= E u ** u fj s « ( 38 ) 

ij 

where the matrix U is given by: 



Qz 


cos 6» 


q z sin^n 


f l li 




sin# |t 


cos 6^ 







cos 9 


q n sin^.. 


Qz 



U = -sin 61, cose,, | (39) 
\q {l cos0n q {{ sin I 

where q z = fa, g y = ^ and <f = , 6*,, ). 
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By using equations l|37fl . H38|) an d E3) we obtain expressions for the matrix 
elements (q,q z ) for different phonon modes A: 



fi|q| 



2pVv L 



fi|q| 



2pVv^ 



rTA2 

Z ±i(±|) 



± 



fi|q| 



(40) 



2pVv T 

where a! and 6' are obtained using a, 6 and c (see equation (|37fl ) and the + and — 
signs on the right hand side of equation (|4(J|) correspond to the light- and heavy-holes 
respectively. A, (i — 1,2,3) are explicit functions of q , q z and 6^, and v TA1 and 
v TA2 are the respective sound velocities corresponding to polarization in the x and 
y-direction. 

The Hamiltonian H^^ ph describing the exciton-acoustic phonon scattering due 
to the piezoelectric effect for LA and TA phonons can be obtained using the Bir-Pikus 
Hamiltonian in equation (|36(l and the matrix U in equation (|39J) . in the form: 



rVPicz 
ex— ph 



M x (q, q z ) \a\ + Sr z ) a n (r' z ) - c\^{a z ) c^a'J 



E 

k,q,qz,r' z ,a> z 

x [b{(-q,-qz)+bx(q,q z ) 
where the matrix elements Mp(q,q z ) are given by: 



(41) 



M^ A (q 7 q z ) =i. 
M^ A1 (q,q z )=i, 
M$ A2 (q,q z )=i, 



2p^Vv, A H 



[3ehuq*q z sin 20^ 



sin 26, 



2 P ^v\ A2 H ( 2e/l14 ^ (cos 



(42) 



where h\± is the piezoelectric coupling constant. Unlike the matrix elements for 
the deformation potential M£ (q,q z ), which depend on the hole spin, the elements 
Mp(q,q z ) corresponding to the piezoelectric interaction are independent of the hole 
spin. 



5. Exciton decay and scattering involving acoustic phonons 

In this section we derive rates of spin transitions due, firstly, to the decay of heavy- 
and light-hole excitons, secondly, to intraband scattering (heavy-hole to heavy-hole, 
light-hole to light-hole) and, thirdly, to interband scattering (heavy-hole to light-hole, 
light-hole to heavy-hole) due to interactions with LA and TA phonons. For the decay 
process, we suppose that initially an exciton with wavevector K and spin S decays 
into an electron hole pair by absorbing an acoustic phonon. We assume that the 
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wavefunction of the final state can be written in terms of a free electron-hole pair as: 

\k, 5) dccay = k-R,n 4^) 4^) b ^ ««) i n «*)' ( 43 ) 

where we have included the allowed A-modes within the summation, although the 
contribution of each A-mode phonon interaction to the decay process can be analysed 
separately by excluding A from the summation, and where n^ qz denotes the occupation 
number of phonons. 

For the scattering process, we suppose that an exciton is scattered from its initial 
state of wavevector K and spin S to a final state with wavevector K' K) and spin 
S' through the following channels which may occur with or without a change of the 
exciton spin: 

(i) \K;S) — ► \K'\S') with energy h\o^_g, A (acoustic phonon emission), 

(ii) \K; S) — > \K' ; S') with energy hujg,_^ A (acoustic phonon absorption). 
The spin relaxation time T sp is calculated using the Fermi golden rule: 



1 



' sp 



2vr 



E E 

k,q,q z 



K/|#mtK>| 2 (h+g + !) hN q x 5(E C x(k + q) - E ex (k) - hu qX ) 



+ {f%_* + 1) f% m + 1) S(E(k -q)- E cx (k) + hLJ qX ) 



(44) 



where \i) and |/) denote the initial and final states, Hi nt is the interaction operator 
and N q is the thermalised average number of phonons at the lattice temperature T\ at , 
and is given by the Bose-Einstein distribution 



N x = 
q 



exp 



- 1 



The exciton distribution function f? is given by 0: 

E CK (k) - [i c 



h 1 = exp 



knT„. 



/csTcx In 



1 — exp 



- 1 



-2irh 2 n ex a e 



(45) 



g m c k B T c: 

where T ex is the exciton temperature, fee is the Boltzmann constant, \i c is the chemical 
potential, n CK is the exciton density, a c is defined in equation 1)25(1 . g is the spin 
degeneracy factor and, as before, m is the effective electron mass. The summation 
in equation (|44|l includes all the initial and final spin states of the coupled electron- 
hole system and thus the relaxation time r sp incorporates the effect of changes from 
spin-up (spin-down) to spin-down (spin-up) electron and hole states of the exciton. 

By using the initial state given by equation ()26|1 . the final state given by 
equation (|43|1 . and the Hamiltonian operator of exciton-phonon interactions due to 
the deformation potential in equation (|32[) . we obtain the transition matrix elements 
for the decay of a heavy-hole exciton into a free electron heavy-hole pair as: 



ph 



(46) 



Af T A , (q, q z ) F c (q z ) V ls {a h K + k h ) - M± 3 _ (q, q z ) F h (q z ) * ls (fc c - a e K) 

d^l^M «* A _ to ^(f^+l)f~ k (N q +S ±hl ) 
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where the signs + and — in <5±i i correspond to a phonon emission and absorption 
process respectively and vice versa in ■ An expression similar to that in equation 
(|4ti|l applies for the decay of the light-hole exciton into a free electron- light-hole pair 



when the reduced rotation matrix d 3 ^ , 3 ((f)) is replaced by d 0/ i , i ((f)) and when the 

fJ-2 M 2 ^2 A< 2 

matrix element M+Aq, q z ) is replaced by M\(q,q z ). Appropriate changes should 

± 2 ± 2 

also be made to a c and ah (which corresponds to the light-hole in equation (J2SJ)) in 
the exciton wavefunction ^i s (a c kh + ahk e ). 

By using equation l|41|) . we obtain the transition matrix elements for the decay 
of the heavy-hole exciton due to piezoelectric effect as: 



3/2 



(f\H, 



Piez i -vHH _ ».rA, 
ex— phl'/c 



Af£(£ &) *i s (a h i^ + fch) F e (q z ) - *i s (fc c - a c i^) F h {q z ) 



xdK 2 ,,(</>) d 



1/2 



An expression similar to this applies for the decay of the light-hole exciton when the 
reduced rotation matrix d 3 ^ , 3 (6) is replaced by d 3 ^ 2 , 1 (</>). 

The transition matrix elements corresponding to the scattering of a heavy-hole 
exciton to another heavy-hole state due to deformation potential interaction is obtained 
as: 

(/l^ex-phK)scatt-DP = E 

Ma 



[-M* (&<fe) F e (g») F c (q' z ) V ls (a h K - k) ^ ls (a h K' - k) 



M^s (q, q z ) F h (q z ) F h (q' z ) * ls (a h K' + k ± q) <£ ls (a h K + k) 

2 

,1/2 



X^| Mi W<f ^W^^ ^/(/^+l)/ E (^+« ± i,i) (47) 



where the reduced rotation matrix d 



"' 2 , 3 ((/)) is replaced by d 3 ^ 2 , i ((f)) and the matrix 

2 ' M 2 ^ 2 



element M^3(q,q z ) replaced by x (g, g z ) for the case of scattering of a light-hole 

± 2 ± 2 

exciton to another light-hole state. For the scattering of a heavy-hole exciton to a 



light-hole state and vice-versa, d 3 ^ 3 2 , 3 ((f>) should be replaced by d a/ 3 , t ((f)) instead. 

^2^2 ^2^2 

The transition matrix elements corresponding to piezoelectric effect are not shown 
here but can be derived in similar manner as equation (|47(l . 

The spin relaxation rates for various phonon mediated processes (such as those 
due to the deformation potential and piezoelectric interactions), are finally calculated 
using equations (|46|) - (|47(l in (|44|l and by using the following summation over spin 
states (see Brink [35], page 24): 



,3/2 



E 

E | 



UK 2 (1 w 



= 2 



j3/2 



((f) =\ (1 + 3 cos 2 



I sin 2 r 



(48) 



6. Numerical Results and Discussion 



In figure [5] we plot the spin relaxation time r sp of the heavy- and light-hole excitons 
as functions of the exciton wavevector \K\ in GaAs7Alo.3Gao.7As quantum wells of 
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well width 70 A and lattice temperature T\ &t — 10 K. We calculate r sp corresponding 
to intraband scatterings of heavy-hole to heavy-hole (labelled HH — » HH) due to the 
deformation potential interaction using equations IMl't , H35[) , (|40|) , (|44|l , 145( 1 and (|47|l 
and integrating over both the in-plane (q) and perpendicular (q z ) components of the 
phonon wavevector. A similar approach is used to calculate r sp due to light-hole 
to light-hole scattering (LH — > LH) via the deformation potential interaction. r sp 
due to the decay of the heavy-hole exciton (HH — > e + h) and decay of light-hole 
exciton (LH — > e + h) is likewise calculated using equations J221>, JUSJ: 
and 1 (45( 1 and the appropriate transition matrix elements. The material parameters 
used in the calculation (following Hj,[II]-|2I]) are listed in tabled We employ the 
BcnDanicl-Duke model JjJlJ and the approach outlined in the Appendix to calculate the 
electron and hole wavefunctions £,g <(z e ) and .(zh) respectively. We have assumed 

v tai ~ v ta2 an d have neglected the effective mass and dielectric mismatches between 
the well and barrier materials in our calculations. 

Our results show that as the wavevector K increases, the exciton spin relaxation 
is likely to occur via the decay channel rather than by the scattering mechanism. 
Depending on K, the heavy-hole and light-hole excitons assume different roles under 
the scattering and decay processes. The light-hole exciton is more likely to decay 
(at large K) while the heavy-hole exciton is more likely to be scattered (at small 
K) during spin relaxation. It is interesting to note that the estimated times for r sp 
are reasonably close for the decay and scattering mechanisms in the range 3 x 10 7 
m _1 ^ A' ^ 5 x 10 7 m _1 . This indicates the delicate balance between decay and 
scattering processes which gives rise to spin relaxation mechanisms and the need for 
a careful examination of the effect of the exciton temperature T cx and the exciton 
density n cx on exciton spin dynamics in low dimensional systems. 

While LA and TA phonons are both effective in bringing about spin relaxation, the 
contribution from TA phonons appears to dominate the spin dynamics of excitons for 
a wide range of K. The spin relaxation times (via the scattering and decay processes) 
are calculated to be of the order of 10 to 100 ps which is consistent with experimental 
results Elj measuring r sp in the order of a 100 picoseconds. In table El we have 
compared r sp (ps) corresponding to various scattering processes in GaAs/Alo.3Gao.7As 
quantum wells at well width 100 A, K = 2 x 10 7 m" 1 and T iat = 4.2 K. The 
calculated values show that spin relaxation is dominated by TA phonons coupled via 
the piezoelectric interaction to excitons undergoing intraband (HH ^ LH) scattering. 

In figure|21we plot the spin transition time r sp as a function of dimensionality a of 
the heavy-hole exciton in GaAs/Alo.3Gao.7As material systems. The contribution of 
both LA and TA phonons is included in the decay and intraband scattering calculations 
in which only the deformation potential interaction is taken into account. Figure 01 
shows clearly the explicit dependence of r sp on a. It is interesting to note that r sp rises 
as high as 1 x 10 3 ps - 1 x 10 6 ps for a ^ 1.3. Experimental results JSj have shown that 
the spin relaxation time of excitons in quantum disks is much higher (w 1000 ps) than 
in quantum wells. Hence at small dimensionalities there is a substantial increase in 
T sp . In table 01 we have compared r sp due to the deformation potential interaction for 
different configurations of GaAs/Al .3Gao.7As systems. The results show the striking 
effect on a and r sp of changing the radius of both the quantum wire and the quantum 
disk and its thickness. 
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Fig. 2 
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Figure 2. The spin transition time r sp as a function of the exciton wavevector 
K in GaAs/Alo.3Gao.7As quantum wells. r sp involves intraband scatterings of 
heavy-hole to heavy-hole (labelled as HH — > HH) and light-hole to light-hole 
scattering (LH — > LH), and decay of the heavy-hole exciton (HH — > e + h) 
and light- hole exciton (LH — > c + h). The temperatures are Ti at = 10 K and 
T ox = 20 K. 
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Fig. 3 
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Figure 3. Spin transition time r sp as a function of the dimensionality a of the 
heavy-hole exciton in the GaAs/Alo.3Gao.7As material system at K = 1 xl0 7 m -1 
and Tj at = 4.2 K, T ex = 20 K. All other parameters used are listed in tabled The 
contribution of both LA and TA phonons are included in the decay and intraband 
scattering calculations. 
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Table 1. Material parameters. 



p 


5.318 x 10 3 kg/m 3 


e 


12.9 


71 


6.85 


V TA 


1.86 x 10 3 m/s 


u 


0.31 eV 


72 


2.1 




4.25 X 10 3 m/s 


b 


2.86 eV 


73 


2.9 


hiA 


0.16 C/m 2 


c 


-2.27 eV 


9 


4 


»— 'C 


-8.0 eV 


T 

± ox 


20 K 




5 x 10 9 



Table 2. Comparison of r sp (ps) corresponding to various scattering processes in 
GaAs/Alo.3Gao.7As quantum well at well width 100 A, \K\ = 2 xl0 7 m _1 , with 
T lat = 4.2 K, T ox = 20 K. 
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Table 3. Comparison of r sp of the heavy-hole exciton due to the deformation 
potential interaction for different configurations of GaAs/Alo.3Gao.7As systems 
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= 1 xlO 7 ra- 
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Radius 


Thickness 


Dimensionality 


Tsp(ps) 


r ap (ps) 




R(A) 


L(A) 


a 


HH^HH (Scat) 


HH->e+h (Decay) 


quantum wire 


50 




1.31 


450 


1 x 10 3 


quantum wire 


100 




1.93 


82 


150 


quantum disk 


30 


50 


0.85 


2.7 x 10 3 


1.1 x 10 4 


quantum disk 


100 


100 


1.72 


135 


223 


quantum disk 


300 


150 


2.24 


50 


8.3 



7. Conclusions 

We have performed a comprehensive investigation of exciton spin relaxation involving 
LA and TA acoustic phonons for various configurations of the GaAs/Al .3Gao.7As 
material system. Our results show the explicit dependence of r sp on a which is an 
important result as the relative motion of an exciton is possibly best described by an 
intermediate dimension in quantum wells, wires or disks. Our results also highlight 
the specific role played by TA phonons coupled via the piezoelectric interaction to 
excitons. Our calculations show good agreement with experimental measurements 
and can be easily extended to coupled double quantum wells as well as to indirect 
excitons and impurities confined in low dimensional systems. Our results show that 
depending on the exciton centre-of-mass wavevector, the balance between decay and 
scattering mechanisms resulting in exciton spin relaxation can be delicate. Future 
developments in this field of research will therefore require a careful examination of 
the conditions (for example the exciton temperature, density and configuration) before 
an accurate picture can be drawn of the spin dynamics of excitons in semiconductor 
systems. 
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8. Appendix 

The hole wavefunction can be written as a superposition of four states of the valence 
band top (labelled bym = ±§,±|), characterized by the Luttinger parameters 71, 72 
and 73 which appear in the 4x4 Luttinger Hamiltonian |28| : 



He = 



c* 
b* 




(49) 



\ 

-b 

c 

Hhh J 

where I4 is the 4x4 identity matrix and where the heavy-hole and light-hole 
Hamiltonians i?hh, -ffm an d the mixing parameters are given by 

. . h 2 (k 2 x + k 2 ) 



c 
H lh 


-b* 
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^lh 

C* 



H hh = - 



Hn 



2mo 
1 



2mo 

h 2 

c(k x , k y ) = VSj 2 7; 

ZTUq 

b(k x ,k y ,p Zh ) = V3j2 



«}(71-272) 
(pi) (71 + 2 72 ) - 



2mo 
t?{kl + 



2mo 



" (7i 

(7i " 



■72) 



72 j 



2mo 



(PzhX^a: -ifcy), 



where m is the free-electron mass and V^z^, L w ) is the square well potential 
of thickness L w that confines the holes, kh = (k Xh , k Vh , k Zh ) denotes the three- 
dimensional wavevector of the holes. We have also used the axial approximation 
72 ~ 73 to simplify c(fch) to a single term and have replaced k Zh by —id/dzh in the 
expression for (j> Zh )- The material parameters 71,72,73 are related to the inverse 
effective masses of electrons and holes. In the k Zh direction the heavy-hole, light-hole 
masses are given by, respectively, 

1 1 

TOhh = 7. — i TO ih = — ■ 

7i - 272 71 + 272 

The transverse masses of the charge carriers are given by the coefficients of the k 2 x + k 2 
terms. Consequently the light-hole transverse mass is larger than the heavy-hole 
transverse mass. 

The hole wavefunctions in a quantum well are obtained by solving: 



E 



(<Jz\H c \cr z ) + Vh(zh, £w)<5o% o-^ Cfc'^h) = Eh(kh) (zh) 



(50) 



where the hole energy E^k^) corresponds to a heavy-hole or light-hole, depending on 
the mass used to model the profile of the potential function Vh(^h)- The summation is 
over the four spin projections a z = ±|,±i of the hole spin, with a aligned along the 



quantization axis. When b 



the Hamiltonian H c splits into two independent 



Hamiltonians, one each for the heavy- hole (a z = ±1) and light-hole (a z = ±5) cases. 
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